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We study the mapping between time-dependent densities and potentials for noninteracting elec- 
tronic systems on lattices. As discovered recently by Baer [J. Chem. Phys. 128, 044103 (2008)], 
there exist well-behaved time-dependent density functions on lattices which cannot be associated 
with any real time-dependent potential. This breakdown of time-dependent U-representability can 
be tracked down to problems with the continuity equation which arise from discretization of the 
kinetic-energy operator. Examples are given for lattices with two points and with TV points, and 
implications for practical numerical applications of time-dependent density-functional theory are 
discussed. In the continuum limit, time-dependent noninteracting U-representability is restored. 



I. INTRODUCTION 

Time-dependent density-functional theory (TDDFT) 
[H is a widely used methodology for the description of 
electron dynamics in atoms, molecules and solids. The 
birth of TDDFT dates back to the seminal work by 
Runge and Gross (RG) @ who proved a one-to-one cor- 
respondence between time-dependent densities, n(r, t), 
and time-dependent potentials U(r, t). The RG theorem 
Q can be viewed as the time-dependent analog of the 
Hohenberg-Kohn theorem Q in static DFT Q. 

While TDDFT has found widespread practical use for 
describing such diverse phenomena as molecular excita- 
tions, nanoscale electron transport, or strong- field pro- 
cesses, there is a continuing need to study the funda- 
mental underpinnings of the theory. Some of the basic 
issues in TDDFT that have recently received much atten- 
tion are nonadiabaticity and memory effects [H, Q , and 
causality and the time-dependent variation principle 0]. 

This paper has been motivated by an unexpected and 
disquieting discovery made recently by Baer [8j: on lat- 
tice systems there exist seemingly well-behaved time- 
dependent density functions which are not Unrepresen- 
table (VR), i.e., which cannot be associated with any 
time-dependent potential. Attempts to invert the time- 
dependent Schrodinger equation to explicitly construct 
the potential from such densities encounter severe numer- 
ical instabilities. It was suggested in Ref. [1] that such 
instabilities might be potentially disastrous for practi- 
cal applications of TDDFT, since they could cause the 
exchange-correlation (XC) potential U xc to become an 
extremely sensitive functional of the density. 

The purpose of this paper is to shed light on two im- 
portant questions which are raised by Baer's study: (1) 
what is the fundamental reason for the absence of time- 
dependent U-rcprcsentability on lattices, and (2) how 
does time-dependent U-representability emerge in the 
continuum limit? Both questions have a profound im- 
pact on our fundamental understanding of TDDFT as 
well as on practical applications. We will give some gen- 
eral answers, as well as some specific illustrations for the 
mapping between time-dependent potentials and nonin- 
teracting densities on one-dimensional lattices. 



II. DFT ON LATTICE SPACES 

The basic existence theorem of ground-state DFT, the 
(first) HK theorem 0], states that there is a 1:1 corre- 
spondence between external potentials and ground-state 
densities. In other words, it cannot happen that two po- 
tentials that differ by more than a constant produce the 
same ground-state density of an Af-particle system. The 
proof of this theorem is based on the Ritz variational 
principle: 

(* gs |ff|* gs ) < (*|#|*>, (1) 

where H is the Hamiltonian of the interacting many- 
particle system, ^/ gs is the associated ground-state wave 
function, and ^Sf is any other many-particle wave func- 
tion. Notice that this fundamental variational principle 
does not require that the system lives in a continuous 
coordinate space; it remains valid for systems that are 
defined on a discrete space, such as on a real-space lat- 
tice or using a finite set of basis functions. 

The basic DFT framework assumes interacting and 
noninteracting U-representability Q. In other words, it 
is assumed that any density function which is mathemat- 
ically reasonable (i.e., does not diverge, and integrates to 
N) is a ground-state density belonging to some external 
potential. This property is of obvious importance for a 
mathematically meaningful implementation of variations 
with respect to the density. Unfortunately, attempts to 
prove U-representability for all densities in the contin- 
uum case (e.g. in Ref. [lOj]) have so far not been fully 
successful. 

In ground-state DFT, lattice systems have been much 
used to address the U-representability problem. In 1983, 
Kohn 

El 

showed that in the vector space of density 
functions that are defined on discrete lattices, a den- 
sity in the neighborhood of a VR density is also VR. 
In 1985, Chayes, Chayes and Ruskai [l^] proved that any 
mathematically well-behaved density function on finite 
or infinite lattices can be represented as the density of 
a pure ground state or as an ensemble density associ- 
ated with a degenerate ground state. In density vector 
spaces of arbitrary dimension, it is equally likely to en- 
counter pure-state and ensemble-state VR densities [l3[ . 



2 



Since the lattice can be arbitrarily dense, this effectively 
solves the V-representability problem in DFT. In a re- 
lated approach, U-representability was recently proved 
by Lammert for coarse-grained systems [13] • 

We mention that lattice systems have also been very 
helpful in studying the issue of nonuniqueness in spin- 
DFT, where it was discovered that there are sev- 
eral different types of situations where the mapping be- 
tween densities and spin magnetizations, {n(r), m(r)}, 
and potentials and magnetic fields, {U(r),B(r)}, is not 
unique. In the continuum limit, only some special cases 
of nonuniqueness in spin-DFT survive. 



III. V-REPRESENTABILITY IN CONTINUOUS- 
AND DISCRETE-SPACE TDDFT 



Let us now treat the time variable as continuous but the 
space variable x as discretized on a grid with equidistant 
grid spacing a. The position of the jth lattice point is Xj, 
and rij denotes the number of particles in a bin of width a, 
centered at Xj. Thus, the lattice quantity corresponding 
to the particle density n(x 7 t) is rij(t)/a. The advantage 
of this definition is that the normalization condition on 
an TV-point lattice becomes simply 



(4) 



For the wave functions, we have ip(x,t) — ► ipj(t)/^/a. 

Using a three-point finite-difference representation of 
the Laplacian operator (which is accurate to within terms 
of order a 2 ) [l7|, we have 



The issue of U-representability in TDDFT turns out 
to be quite different from static DFT. Van Leeuwen (vL) 
[l6| showed that, assuming that the initial state is non- 
interacting VR, there is always a unique time-dependent 
potential V(r, t) in a noninteracting system that pro- 
duces a given interacting density n(r, t) at all times. 
In other words, the time-dependent U-representability 
problem in the continuum case can be considered solved 
by the vL-construction. In view of this, the find- 
ing by Baer @ that time-dependent noninteracting V- 
representability breaks down on lattices is particularly 
unsettling and calls for an explanation. 

Unlike the HK theorem in static DFT, the RG theo- 
rem and the vL construction are not based on 
a minimum principle. Instead, it is demonstrated that 
potentials which differ from each other by more than 
a time-dependent function cause the time evolution of 
the respective systems to proceed in such a manner that 
the time-dependent densities are not the same. The 
original RG proof requires two steps: first it is shown 
that if a system evolves from a given initial state un- 
der the influence of two external potentials V(r, t) and 
V'(r,t) 7^ V(r,t) + c(t) then the current densities j(r, t) 
and j'(r, t) are different. The next step shows by virtue 
of the continuity equation that the associated particle 
densities n(r,t) and n'(r,t) become different, too. The 
continuity equation also plays a crucial role in the vL 
construction. 

Let us therefore take a closer look at the continuity 
equation, given by 



n(r,t) = -V-j(r,t) ! 



(2) 



where the dot denotes a partial derivative with respect to 
time. Here and in the following, we use Hartree atomic 
units with K = e = m = 1. For the simplest case, a 
single particle in one dimension, the continuity equation 
becomes 



n(x, t) = 



r{x,t) 



d 2 ?/>(x,t) 
dx 2 



ip{x,t) 



dip*(x,t) 
dx 2 



(3) 



i 

2a 2 



(5) 



where the time arguments are suppressed for brevity. 
Without loss of generality we take the following form of 
the wave function: 



rij e 



(6) 



where aj is a real-valued time-dependent phase. The 
discretized continuity equation then becomes 



1 



hj = -^y/n 3 n j+1 sin( aj - a j+ i) 
1 , . , 



(7) 



We will now show that hj is not allowed to take on arbi- 
trary values. We have 



1 



\hj\ < -J [y/n jnj+1 + y/rijUj. 



(8) 



Next, let nj + i = rij + A. Using the normalization condi- 
tion ([!]) one obtains rij—\ < 1 — 2rij — A (the equal sign 
arises for lattices consisting only of three points, or if the 
potential is infinity on all but the three points j — 1, j, 
and j + 1). Thus, 



N < 



^rij + A + y/l - 2n~j 



(9) 



The right-hand side is maximized for A = and rij = 
(3 + V3)/12, which leads to 



M < 



0.683 



(10) 



On a one-dimensional, equidistant grid with a three-point 
finite-difference representation of the kinetic energy op- 
erator, the time derivative of the density is restricted by 
the local upper bound (fTU]) . In other words, the density 
on a given grid point cannot change arbitrarily fast. 
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The upper bound (fTU)) can be generalized for a d- 
dimensional grid using a fc-point formula for the Lapla- 
cian operator [13], resulting in \ fij\ < Af/a 2 , where Af is 
a finite numerical constant. It is easy to see that on an 
equidistant <i-dimensional grid, Af. — dA\ . Going beyond 
three-point formulas, one finds for instance Af = (4/3)^4|f 
for five-point formulas, and A^ — > 2A% for large k. 

Thus, there are fundamental limits, determined by the 
grid spacing, the number of dimensions, and the discrete 
representation of the kinetic-energy operator, on how fast 
lattice densities can locally change without violating the 
continuity equation. The RG theorem and the vL con- 
struction are thus restricted to time-dependent densities 
which satisfy these bounds. Otherwise, time-dependent 
1^-representability cannot be guaranteed. We will now 
illustrate this with some explicit examples. 



IV. THE 2-POINT LATTICE 

The case of the 2-point lattice was already treated by 
Baer Q; nevertheless, it is a very instructive and peda- 
gogical model, and we will discuss it here in some detail. 

We start with the single-electron time-dependent 
Schrodinger equation (TDSE) in real space: 



«^(r,t) 



+ V(r,t) 



M,t) 



(11) 



On a two-point lattice, the wave function is completely 
specified by its values on point 1 and on point 2, ipi(t) and 
tp 2 {t). These two (complex) quantities can be arranged 
as the components of a vector. The TDSE can then be 
discretized and written in the form of a 2 x 2 matrix 
equation: 



d 



Mt) 
Mt) 



i 

"To 1 




where we used a standard 3-point finite-difference rep- 
resentation of the kinetic-energy operator assuming that 
the two lattice points have a spacing a. 

The notation can be simplified recognizing that an 
overall constant shift of the potential is irrelevant. This 
means that the potential can be redefined as Vj = -^r+Vj, 
j = 1,2, and 



d 



Mt) 



Vi(t) 



Mt) 
Mt) 



(13) 



We now write the wave function on points 1 and 2 in the 
form ©. The densities and phases, ni, 2 and 0^2, are 
real functions of time. Inserting this in the TDSE (fT3"|) 
gives the following equations for the real and imaginary 
part of the potential on point 1 : 



$tVi = -ai + TTTA/— cos(a 2 -ai) 
2a z V n\ 



hi 1 pna . , 
2^ + 2^Vm Sm(a2 



(14) 
(15) 



and similar on point 2. If we require the potential to be 
purely real, then the imaginary part, given by equation 
(fl~5|) . vanishes. We can then solve for the phase difference: 



OL\ — Q.1 



a ni 



(16) 



Due to norm conservation ((4]) , we have only one indepen- 
dent density variable, since n\ — 1 — n 2 . Likewise, there 
is only one independent potential variable, since an ar- 
bitrary overall constant shift only affects the phases 0]. 
Thus, only the potential difference matters, and we let 
AV = nV 2 - SRVi. From Eq. we get 



AV 



(ai2 - ai J 



1 

2^2 




4/ — cos(ai — q.%) 



and the phase differences can be eliminated using Eq. 
(fT5)) . This yields the potential difference as a function of 
the densities at point 1 and 2: 



Ay = 



\Jn\n-2 - 
ni - n 2 
2a 2 niri2 



"1 



h\{n\ - na) 



(18) 



This equation for AV is well-behaved and has a unique 
solution as long as niri2 — a A h\ > 0, or 



K,2 < 2 



(19) 



[n\ and n 2 themselves are of course always > 0; further- 
more, hi = —hi)- This yields the following upper bound: 



0.354 



"1.2 



(20) 



Time-dependent single-particle lattice densities Ui,ri2 
which violate conditions (fT9|) or (|20|) at any time t do 
not produce a real value of AV at that time and are 
therefore not time-dependent VR. This upper bound for 
the two-point lattice is clearly more restrictive than con- 
dition (jTUJ) for a general one-dimensional lattice. 

As an illustration, let us consider a one-dimensional 
time-dependent density function of the form 

n(x,t) = A\ cos 2 (7rx) + A\ sin 2 (27ra) 

+ 2AiA 2 cos(7rx) sin(27ra) cos(ujt) , (21) 

— 1/2 < x < 1/2, arising from a superposition of the 
first and second single-particle eigenstate of a particle 
in a one-dimensional box of length one. We discretize 
n(x,t) on a 2-point lattice such that ni j2 = n(xi_ 2 ,t), 
where Xi — —1/6 and x 2 = 1/6, and the lattice spacing is 
a = 1/3. The lattice normalization condition Q requires 
choosing Ai and A 2 such that rii + n 2 = 1. 

The top panel of Fig. [1] shows ni and n 2 as a func- 
tion of time for a frequency uj = O.Itt and A\ = 0.77 
and A 2 = 0.27, together with the associated potential 
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10 20 30 40 50 60 70 80 90 100 

t (a.u.) 

FIG. 1: Top: two-point lattice densities n\ (full line) and ri2 
(dashed line) resulting from n(x, t) of Eq. (121[) . with uj = 0.17T. 
Bottom: potential difference AV = Vi - V [Eq. JTHJl]. 



difference Ay, following from Eq. (|T8|) . This is a case 
where the lattice densities are VR, and the potential is 
a well-behaved function of time which drives the particle 
density periodically from the first to the second lattice 
point and back. In fact, in this low-frequency case the 
adiabatic approximation for the potential, 



m(t) - n 2 (t) 



2a? y /n 1 (t)n 2 (t) 



(22) 



(which assumes that at each moment in time the system 
is in the ground state associated with the instantaneous 
potential), is indistinguishable from the exact AV(t). 

Figure [2] shows the potential reproducing ri\ and n 2 
at higher frequencies, from uj = tt through 57r. As u> 
grows, the time-dependent potential AV starts to dif- 
fer more and more from the adiabatic approximation 
AVadia- In fact, there is a crossover at uj — 2.8657T 
where one finds AV(t) = 0, which corresponds to free 
charge-density oscillations on the two-point lattice with 
frequency u> = 1 /a 2 (the difference of the two eigenvalues 
of the static Schrodinger equation with V\ = V> = 0). For 
lu > 2.8657T, AV and AVadia are out of phase, which may 
be viewed as an indication of entering the high-frequency 
regime. 

For uj > 4.5877T one finds that AV diverges period- 
ically, and there are regions in time where there is no 
solution at all, as seen for u> = 5ir in panel (e) of Figure 
[2j This is a striking example where the time-dependent 
lattice density is non-VR. 

In the regions in which AV does not exist, condition 
(fT9)l is violated, which implies that the densities change 
too rapidly. In other words, there is a limit on how 
much density can move how fast between the two points. 
Violation of this limit means that no real-valued time- 
dependent lattice potential can be found which is capable 
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FIG. 2: Full lines: potential difference AV = V 2 - V for 
uj — n, 2n, 3n, 47r, and 5n (a-e), with T = 2n/oj. Dashed lines: 
adiabatic potential [Eq. (122[l ]. Panel (e) shows regions where 
no potential exists that produces the given lattice densities 
ni(t) and ri2(t), which means that they are non-VR. 



of driving the system strongly and fast enough to produce 
the given time-dependent density. 

To conclude this section, we mention that some of the 
issues addressed here have also been recently discussed 
by Verdozzi in the context of Hubbard dimers [l8| . 



V. THE iV-POINT LATTICE 

For a linear TV-point lattice with equidistant grid spac- 
ing a and a 3-point finite-difference formula, the dis- 
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crctized TDSE reads 
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and V; 



(23) 
have been 



where the time arguments of the ip 
omitted and the constant I/a 2 coming from the finite- 
difference kinetic-energy operator has been absorbed in 



the potentials Vi(t), as in Eq. ([II 

It is straightforward to generalize the analysis of the 
previous section to the TV-point lattice. We again use 
ansatz ([6]) and eliminate the phases using the requirement 
that the potential be real at each point. Defining 



(24) 



we obtain after some algebra the following expression for 
the point-to-point potential differences: 



V k+1 - V k 



(n k h k+ i + n k+ ih k )S k 
2n k n k+1 



n k -n k+1 , 2 



+ ^n k+1 n k+2 - a±Sl +1 - ^-Jn k ^n k - a^U , 



,N-1. 



(25) 



For notational convenience, Eq. (|25j) makes reference to 
the zeroth and the (TV + I)st lattice point. These points 
should be regarded as "virtual points" , consistent with 
the boundary condition tiq = n^ + i = 0. With this and 
the normalization condition (0} it is easy to see that Eq. 
(pjj) reduces to Eq. (fT5|) for TV = 2. 

Eq. ([23)) yields real solutions for the point-to-point 
potential differences Vfe+i — V k as long the terms under 
the square roots remain positive. This immediately leads 
to the following constraint on the time-dependent lattice 
density, which is a generalization of Eq. : 

Vn^T^ fe=l,... 5 AT-l. (26) 

This condition represents a much tighter and more spe- 
cific criterion for the V-representability of a given time- 
dependent lattice density than the upper bound (fTU)) . 

Let us compare the two criteria for the same oscillating 
density function (JUI), with A x = 0.77 and A 2 = 0.27, 
that was studied in section HVl We consider TV-point 
lattices with grid spacing a — 1/(N+1) and lattice points 
at Xj = ja — 1/2, j = 1, . . . ,N, and we normalize the 
discretized density using 

The full line in figure [3] shows the critical frequency at 
which the time-dependent density becomes non-VR on at 
least one point on an TV-point lattice, according to Eq. 
(f2"6")l . By comparison, the dashed line shows the frequency 
at which the same time-dependent density would start to 
violate the upper bound (|I0|) . One finds that the former 
increases linearly with the number of lattice points N, 
whereas the latter grows as TV 3 . This drastically different 
behavior is not surprising: the absolute upper bound (If 01) 
corresponds to the extreme limit where the density is 
completely concentrated on three lattice points only. For 
N < 4, this gives a quite good agreement as seen in 
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FIG. 3: Full line: critical frequency at which n(x,t) from 
Eq. (|21[) becomes non-VR on an TV-point lattice. Dashed 
line: frequency at which the density would violate the upper 
bound (|10l) . The dots indicate the cases studied in Figure [4] 

figure [31 but the more points one adds to the lattice, the 
less this extreme scenario applies to the density that we 
actually consider here. 

Figure 2] gives further evidence of the breakdown of 
T^-representability. We plot the time-dependent force as- 
sociated with the potential differences ([25]) . 

F k+ i =-(V k+1 -V k )/a, (27) 

which represents a discretized version of the continuum 
expression F(x) = —dV{x)/dx. At TV = 50 and to = 80w 
the density ([2"T[) is VR, as can be seen from Figure [3] (the 
dot below the full line). The forces F k+ i are therefore 
well behaved, i.e., they exist and are finite everywhere on 
the lattice at all times. If we keep the number of lattice 
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FIG. 4: Snapshots of the time-dependent force (|27[) generating 
the density ([21]) for i = (open circles), T/8 (full circles), T/4 
(diamonds), 3T/8 (full squares), and T/2 (empty squares). 
Top panel: force exists everywhere at N = 50 and uj — 80-7T, 
n(x, t) is lattice- VR. Other panels: V^-representability breaks 
down for increasing u at fixed TV, or decreasing TV at fixed ui. 



points fixed but increase the frequency, we cross the line 
of critical frequencies at u> = 1097r: U-representability 
breaks down and there are regions on the lattice and 
times during the cycle where no force exists. A similar 
breakdown happens for fixed frequency if the number of 
lattice points is reduced. 



VI. DISCUSSION 

From the above analysis of one-electron lattice densi- 
ties, we can draw a number of conclusions about nonin- 
teracting U-representability in TDDFT. 

Continuum limit. - Condition (|26|) . as well as the 
upper limit (|10p , both become more and more easily 
satisfied as the grid spacing a decreases. This is also 
evident from Figure [3] Thus, in the continuum limit 
a — > 0, which corresponds to a lattice whose distribu- 
tion of grid points is infinitely dense, all time-dependent 



single-particle densities become VR. 

There are various ways to see how time-dependent V- 
representability emerges in the continuum limit. First of 
all, as discussed in section UTTl the continuity equation is 
in general not valid on lattices, which breaks crucial links 
in the chains of arguments that constitute the RG proof 
and the vL construction. In the a — > limit, the continu- 
ity equation is restored, and thus the fundamental proofs 
of TDDFT which establish the one-on-one mapping be- 
tween densities and potentials. 

Another way of looking at things is to consider the 
energy eigenvalue spectrum associated with a particular 
lattice. Each excitation energy can be associated with a 
resonant electronic eigenmode, i.e., a free (undriven) os- 
cillation of arbitrary amplitude. If the time dependence 
of a given lattice density is close to that of one of the 
lattice eigenmodes, or a superposition thereof, then it is 
possible to find an external time-dependent driving po- 
tential which gives rise to this density. On the other 
hand, if the given time-dependent density changes much 
more rapidly than any available lattice mode, then no ex- 
ternal potential exists which is capable of inducing these 
rapid changes. This was clearly demonstrated in section 
IIVI for the 2-point lattice. 

Ultimately, it is the time-energy uncertainty relation 
At AE > h whi ch g overns how fast a quantum state can 
change in time [19(. Here AE indicates a spectral mea- 
sure of the initial Hamiltonian of the system, such as 
its "energy spread" . Consider a 1-dimensional Appoint 
lattice with grid spacing a, and uniform (constant) po- 
tential. There are N eigenvalues £„, n — 1,...,N, of 
the static Schrodinger equation, and N ~ 1 eigenmodes 
with frequency u n = e n — E\. The dynamical range or 
energy spread of this lattice is given by the difference 
of the highest and the lowest eigenvalue, which scales as 
I /a 2 . According to the uncertainty relation, the range 
of possible response times of the system decreases at the 
same rate. Thus, the smaller the grid spacing, the more 
rapidly the system can change in time, and the wider the 
range of time-dependent U-representability. 

Complex potentials. - It may be of interest to observe 
that non-VR lattice densities can still be associated with 
external potentials. However, these potentials have to 
be complex, and are therefore non-Hermitian. Complex 
potentials are a well-known tool to generate absorbing 
boundary conditions for quantum dynamics simulations 
on finite grids [20]. However, complex time-dependent 
potentials can also be constructed to be norm-conserving, 
which means that the resulting time evolution is unitary. 
This is explicitly demonstrated in the Appendix. 

The physical mechanism by which a complex potential 
reproduces a non-VR density is not by "driving" the den- 
sity from one lattice point to another, as a real potential 
would try to do. Rather, complex potentials introduce 
local particle sources and sinks. A time-dependent den- 
sity that moves rapidly between two lattice points can 
thus be viewed as being "destroyed" on the first point 
and "re-created" on the second point. There is in princi- 
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pie no limit as to how fast such processes can take place. 
This is also reflected in the fact that the continuity equa- 
tions ([2|) and (J3j) are modified for complex potentials to 
include extra source terms. 

The mapping between real densities (a single variable) 
and complex potentials (two variables: real and imagi- 
nary part) obviously cannot be unique, as shown in the 
2-point lattice example in the Appendix. It may be possi- 
ble to restore uniqueness of the mapping by including the 
phase information of the wave function, or equivalcntly 
by using the current as an additional variable. 

Many-electron systems. - In this paper, we have only 
studied the case of a single electron in detail. However, 
the essential arguments for and against noninteracting 
F-representability on lattices, see section Hill can be car- 
ried over to the case of many electrons. In the time- 
dependent Kohn-Sham (TDKS) scheme, the continuity 
equation holds for each individual Kohn-Sham orbital, 
and it is immediately seen that this can be violated on a 
finite lattice. The upper bound p0[) can be generalized to 
the case of several electrons in a straightforward manner. 
It is thus clear that time-dependent y-representability 
for many-electron system is only guaranteed in continu- 
ous space, and not on lattices. 

It is possible to invert the TDKS equation for two elec- 
trons in a doubly occu pied orbital, which was recently 
done for Hooke's atom [2 lL I22I |23| and two-dimensional 
quantum strips [24J. However, for more than two elec- 
trons the construction of the potential in terms of the 
density can no longer be carried out analytically [such as 
in Eq. I|25p]. but only numerically. 

Ensemble-VR densities. - In static DFT, examples 
of non-VR densities were first constructed by Levy [25j | . 
The conundrum was soon resolved by recognizing that 
such densities can come from ensembles of degener- 
ate states. This finally led to the modern view of V- 
representability in DFT in terms of pure- and ensemble- 
state VR densities, see the discussion in section ILT1 

The role played by degeneracies in TDDFT has so far 
not been studied in much detail; the situation is much 
more complex than in static DFT. However, it is highly 
unlikely that the absence of IZ-representability on lat- 
tices can be cured by the same trick that worked in static 
DFT, that is, by representing non-VR densities via en- 
sembles of degenerate states. The simple reason is that, 
unlike in static DFT, non-VR time-dependent densities 
occur already for one-electron systems. Furthermore, the 
physical origins of non-V-representability in TDDFT are 
fundamentally tied to the dynamics (such as the conti- 
nuity equation and how fast quantum states can change 
- see the above discussion). 

Practical consequences. - In Ref. Q it was argued that 
V xc is likely to be an extremely sensitive functional of 
the time-dependent density, especially in strongly time- 
dependent problems, due to the possible occurrence of 
instabilities and breakdown of the density-on-potcntial 
mapping on lattices. Indeed, V xc is defined as a functional 
of the time-dependent density on the domain of densities 



that are noninteracting VR, and is therefore undefined for 
non-VR densities. However, does this have any practical 
consequences for TDDFT? 

V xc needs to "know" whether a time-dependent lattice 
density is noninteracting VR or not. Continuum-space 
XC functionals are not required to have this property, 
since all time-dependent continuum densities are nonin- 
teracting VR by the vL construction [l6| (as long as the 
initial state is VR). Therefore, the exact or any approx- 
imate V xc from continuum TDDFT, if used in a lattice 
TDKS calculation, remains well defined (real and finite) 
even if evaluated with non-VR lattice densities. To deal 
properly and consistently with non-VR lattice densities, 
Vxc must be lattice dependent, i.e., it needs to be con- 
structed for many-body systems on specific lattices. 

Let us now assume that we have the exact V xc for a 
specific lattice system, and we want to carry out a self- 
consistent TDKS calculation using an iterative procedure 
such as discussed in Ref. (2(|. The first iteration step 
involves evaluating V xc with a time-dependent trial den- 
sity. If this trial density turns out to be non-lattice- VR, 
then the TDKS calculation breaks down. If, on the other 
hand, the trial density is lattice- VR, then each subse- 
quent iteration produces only lattice- VR densities, and 
the TDKS calculation will be well-behaved. 

Problems could still arise if V xc were extremely sensi- 
tive within the VR region: for example, if it fluctuated 
wildly for VR lattice densities that are close to the non- 
VR region. However, the examples which we have studied 
here indicate that this is unlikely to be the case. Figure 
2] shows that the force, and hence the potential itself, re- 
mains finite and well-behaved as long as the density is 
VR, and otherwise simply ceases to exist. 

Finally, it is an obvious statement that in the usual 
applications of TDDFT the density follows from the po- 
tential, not the other way round, and therefore remains 
VR by definition. Furthermore, we have seen that non- 
y-representability becomes increasingly rare in the space 
of time-dependent lattice densities if the grid spacing is 
made smaller and smaller by adding more and more lat- 
tice points. We may therefore conclude that the fact that 
the density-on-potential mapping on lattices can break 
down has no practical consequences for TDDFT. 
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APPENDIX: COMPLEX POTENTIALS 

Let us first revisit the derivation of section HVl for the 
two-point lattice, but now admitting a time-dependent 
potential with nonvanishing imaginary part. We impose 
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norm conservation: 

o = h x + h 2 = + M4) + WSife + M4) • (A.i) 

Making use of the TDSE for and ip 2 , this leads to 
the following condition for the imaginary parts of the 
potentials on points 1 and 2: 

mSVi + n 2 ^V 2 = 0. (A.2) 

From Eq. (Ti"5|) we obtain 

ai -a 2 = S m- 1 (^L-2a\f^^V 1 ). (A.3) 
\y<nrnE \ n 2 J 

Taking the time derivative of this and plugging it into 
equation (JT7J) gives, after some algebra, 



\Jn\n 2 — a 4 (hi — 27ii3Vi) 2 

J . h\(nx-n 2 ) n 2 rn \ 

x < ni H 2~ - 2ni3fVi > 

L 2nin 2 n\ J 

+ IV " 2 ^"1^2 - a 4 (»i - 2n 1 5F 1 ) 2 , (A.4) 

which is a generalization of equation (|18p for nonvanish- 
ing imaginary part of the potential. We saw in section ITVl 
that the density is non-VR if the argument n\n 2 — a h\ 
under the square root in Eq. (fl"8|) becomes negative. If 
we admit complex potentials, we have the freedom to 
choose SJVi in such a way that 

rum - a 4 (ni - 2n 1 $V 1 ) 2 > . (A.5) 



This determines $tV 2 — SRFi and ^SV 2 via equations (|A.2[) 
and (|A.4[) . Clearly, 3?Fi is not unique, but each choice 
for QVi which satisfies condition (jA.5|) will reproduce the 
given time-dependent density via the lattice-TDSE (|13|) . 

The arguments presented here for the 2-point lattice 
can be extended for the Appoint lattice in a straightfor- 
ward manner. Generalizing equation (|A.2[) , we find 



N 

5>^=0 (A.6) 



as a necessary condition that the complex lattice poten- 
tial is norm-conserving. The point-to-point differences of 
the real parts of the potential, diVk+i — JiVfc, can then be 
obtained from Eq. (|25|) by replacing Sk in Eq. (|24l) by 



k 

S k =Y J (n -2n 3 ^V 3 ). (A.7) 



We can therefore always satisfy the constraint (|2l)|) 
through appropriate choices of the imaginary part of the 
lattice potential, consistent with the norm-conservation 
condition (| A.6[) . As for the 2-point lattice, these choices 
are not unique. 
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